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Abstract 

Permeability  of  fabric  preforms  and  its  changes  due  to  various  modes  of  the  fabric  distortion  or  deformation  as  well  due  to  fabric  layers 
shifting  and  compacting  is  one  of  the  key  factors  controlling  infiltration  of  the  preforms  with  resin  within  the  common  polymer-matrix 
composite  liquid-molding  fabrication  processes.  While  direct  measurements  of  the  fabric  permeability  generally  yield  the  most  reliable 
results,  a  large  number  of  the  fabric  architectures  used  and  numerous  deformation  and  layers  rearrangement  modes  necessitates  the 
development  and  the  use  of  computational  models  for  the  prediction  of  preform  permeability.  One  such  model,  the  so-called  lubrication 
model,  is  adapted  in  the  present  work  to  study  the  effect  of  the  mold  walls,  the  compaction  pressure,  the  fabric-tows  shearing  and  the  fabric 
layers  shifting  on  permeability  of  the  preforms  based  on  orthogonal  balanced  plain- weave  fabrics.  The  model  predictions  are  compared  with 
their  respective  experimental  counterparts  available  in  the  literature  and  a  reasonably  good  agreement  is  found  between  the  corresponding 
sets  of  results. 

©  2004  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Over  the  last  two  decades,  processing  of  high-performance 
polymer-matrix  composites  via  the  use  of  modern  resin- 
injection  technologies  has  made  major  advances  and  ex¬ 
panded  from  its  aerospace  roots  to  military  and  diverse 
civil  applications.  At  the  same  time,  processing  science  has 
become  an  integral  part  of  the  composite-manufacturing 
technology  so  that  empiricism  and  semi-empiricism  have 
given  way  to  greater  use  of  computer  modeling  and  sim¬ 
ulations  of  the  fabrication  processes.  Among  the  mod¬ 
ern  polymer-matrix  composite  manufacturing  techniques, 
liquid-molding  processes  such  as  resin  transfer  molding 
(RTM),  vacuum  assisted  resin  transfer  molding  (VARTM) 
and  structural  reaction  injection  molding  (SRIM)  have  a 
prominent  place.  A  detailed  review  of  the  major  liquid 
molding  processes  can  be  found  in  the  recent  work  of  Lee 
[1].  One  common  feature  to  all  these  composite  fabrica¬ 
tion  processes  is  the  use  of  low-pressure  infiltration  of  the 
porous  fabric  preforms  with  a  viscous  fluid  (resin).  Con¬ 
formation  of  the  fabric  preforms  to  the  ridges  and  recesses 
in  the  mold  and  the  applied  pressure  can  induce  significant 
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distortions  and  deformations  in  the  fabric  as  well  as  give 
rise  to  shifting  of  the  individual  fabric  layers  and,  in  turn, 
cause  significant  change  in  local  permeability  of  the  pre¬ 
form.  Since  the  infiltrating  fluid  follows  the  path  of  least 
resistance,  local  changes  in  the  fabric  permeability  can  have 
a  great  influence  on  the  mold  filling  process  influencing 
the  filling  time,  the  filling  completeness  and  the  formation 
of  pores  and  dry  spots.  Hence,  the  knowledge  of  the  effect 
of  various  distortion,  shearing  and  shifting  modes  on  the 
fabric  permeability  is  a  critical  step  toward  better  designs 
and  control  of  the  liquid  molding  processes. 

Permeability  of  a  porous  medium  is  one  of  the  most  im¬ 
portant  parameters  controlling  the  flow  of  a  fluid  through 
such  medium.  In  simple  terms,  permeability  can  be  defined 
as  a  (tensorial)  quantity  which  relates  the  local  velocity  vec¬ 
tor  of  the  fluid  flow  with  the  associated  pressure  gradient. 
In  polymer-matrix  composite  liquid-molding  manufactur¬ 
ing  processes  (e.g.  in  the  RTM  and  the  VARTM  processes), 
the  porous  medium  consists  of  woven-  or  weaved-fabric 
preforms  placed  in  the  mold  and  the  fluid  flow  of  in¬ 
terest  involves  preform  infiltration  with  resin.  Complete 
infiltration  of  the  preform  with  resin  is  critical  for  obtain¬ 
ing  high-integrity,  high-quality  composite  structures.  The 
knowledge  of  the  preform  permeability  and  its  changes 
due  to  fabric  bending,  shearing,  compression,  shifting,  etc. 
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Nomenclature 

/  fiber  volume  fraction 

h  fabric  thickness  (m) 

K  permeability  tensor  of  the  fabric  (m2) 

L  length  of  the  quarter  unit  cell  (m) 

p  pressure  (Pa) 

r  fiber  radius  (m) 

s  relative  shift  of  the  adjacent  fabric  layers  (m) 
u  v-component  of  the  resin  velocity  (ms-1) 

U  in-plane  resin  velocity  magnitude  (ms-1) 
v  y-component  of  the  resin  velocity  (ms-1) 

w  z-component  of  the  resin  velocity  (ms-1) 

W  transverse  resin  velocity  magnitude  (ms-1) 

Greek  symbols 

0  resin  viscosity  (N  s  m-2) 

0  shear  angle  (°) 

0  relative  dimensionless  shift  of  the  adjacent 
fabric  layers 

Subscripts 

bot  quantity  associated  with  the  bottom  surface 
of  the  fabric 

corr  quantity  corrected  for  the  effect  of  shear  on  the 
fiber  volume  fraction 

low  quantity  associated  with  the  lower  mold 
surface 

0  quantity  associated  with  un-sheared  fabric 
preform 

top  quantity  associated  with  the  top  surface  of  the 
fabric 

upp  quantity  associated  with  the  upper  mold 
surface 

6  quantity  associated  with  sheared  fabric 
preform 

Superscripts 

B  quantity  associated  with  the  bottom  channel 

F  quantity  associated  with  the  fabric 

T  quantity  associated  with  the  top  channel 


is  crucial  in  the  design  of  a  composite  fabrication  process 
(e.g.  in  the  design  of  the  tool  plate,  or  for  placement  of  the 
resin  injection  ports).  In  general,  the  most  accurate  value 
of  permeability  of  a  porous  medium  is  obtained  by  direct 
experimental  measurements.  However,  the  number  of  fab¬ 
ric  architectures  can  be  quite  large  and  fabric  distortion 
modes  numerous  making  permeability  determinations  via 
the  purely  experimental  means  not  a  very  appealing  alterna¬ 
tive.  In  addition,  sometimes  the  experimentally  determined 
permeability  values  reported  by  different  researchers  for  the 
apparently  identical  fabric  architectures  can  differ  signif¬ 
icantly.  Consequently,  significant  effort  has  been  invested 


Fig.  1.  A  schematic  of  (a)  the  top  view  and  (b)  the  edge  view  of  a 
one-layer  orthogonal  plain-weave  fabric  preform. 


over  the  past  decade  to  develop  computational  tools  which 
can  be  used,  in  conjunction  with  experimental  measure¬ 
ments,  to  determine  the  preform  permeability. 

For  the  computational  modeling  approach  to  be  successful 
in  predicting  permeability  of  the  fabric  preforms,  it  must 
include,  in  a  correct  way,  both  the  actual  architecture  of 
the  fabric  and  the  basic  physics  of  the  flow  through  it.  A 
schematic  of  the  relatively  simple  orthogonal  plain-weave 
one-layer  fabric  architecture  is  shown  in  Fig.  1.  As  seen  in 
Fig.  1,  the  fabric  consists  of  orthogonal  (warp  and  weft)  fiber 
yarns,  which  are  woven  together  to  form  an  interconnected 
network.  Each  yam,  on  the  other  hand,  represents  a  bundle 
of  the  individual  fibers  held  together  with  thread.  In  addition, 
the  fabric  involves  a  network  of  empty  pores  and  channels. 
When  a  fabric  like  the  one  shown  in  Fig.  1  is  being  infiltrated, 
the  resin  flows  mainly  through  the  pores  and  the  channels. 
However,  since  the  fabric  tows  are  porous  (pores  and  channel 
on  a  finer  length  scale  exist  between  the  fibers  in  tows),  the 
resin  also  flows  within  the  yarn.  Thus  when  predicting  the 
effective  permeability  of  a  fabric,  the  computational  model 
must  account  for  both  components  of  the  resin  flow. 

Prediction  of  the  permeability  of  porous  medium  has  been 
the  subject  of  intense  research  for  at  least  last  two  decades. 
Due  to  space  limitations  in  this  paper,  it  is  not  possible  to  dis¬ 
cuss  all  the  models  proposed  over  this  period  of  time.  Nev¬ 
ertheless,  one  can  attempt  to  classify  the  models.  One  such 
classification  involves  the  following  main  types  of  models 
for  permeability  prediction  in  the  porous  media:  (a)  the  phe¬ 
nomenological  models  based  on  the  use  of  well  established 
physical  concepts  such  as  the  capillary  flow,  e.g.  [2,3],  or  the 
lubrication  flow,  e.g.  [4].  These  models  generally  perform 
well  within  isotropic  porous  media  with  a  simple  architec¬ 
ture.  (b)  The  numerical  models  which  are  based  on  numer¬ 
ical  solutions  of  the  governing  differential  equations.  These 
models  generally  attempt  to  realistically  represent  the  archi¬ 
tecture  of  the  fiber  preform  but,  due  to  limitations  in  the 
computer  speed  and  the  memory  size,  are  ultimately  forced 
to  the  introduction  of  a  number  of  major  simplifications,  e.g. 
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[5,6].  (c)  The  models  which  are  based  on  a  balance  of  the 
fabric  architecture  and  the  flow  physics  simplifications,  en¬ 
abling  physically  based  predictions  of  the  preform  perme¬ 
ability  within  reasonably  realistic  fabric  architectures.  One 
of  such  models  is  the  one  proposed  by  Simacek  and  Advani 
[7] .  The  model  of  Simacek  and  Advani  [7]  also  includes  the 
effect  of  important  factors  such  as:  (a)  the  flow  within  the 
fiber  yarn;  (b)  nesting  in  multi-layer  fabric;  and  (c)  distor¬ 
tion  and  deformation  of  the  fabric.  In  the  present  work,  the 
S  model  of  Simacek  and  Advani  [7]  is  extended  to  include 
the  effect  of  shear  of  the  fabric  tows  on  the  effective  volume 
fraction  of  fibers. 

The  organization  of  the  paper  is  as  follows:  a  brief 
overview  of  the  model  proposed  by  Simacek  and  Advani 
[7]  and  its  modifications  are  presented  in  Section  2.  The  ap¬ 
plication  of  their  model  to  reveal  the  role  of  various  fabric 
distortion  and  layers  compaction  phenomena  is  presented 
and  discussed  in  Section  3.  The  main  conclusions  resulted 
from  the  present  work  are  summarized  in  Section  4. 


2.  Computational  procedure 


(width-to-height)  aspect  ratio  (5  or  larger).  The  geometry  of 
the  tows  within  the  cell  can  be  described  using  various  math¬ 
ematical  expressions,  e.g.  [8],  for  the  top,  ztop(X  y),  and  the 
bottom,  Zbotfe  y),  surfaces  of  the  fabric,  respectively.  In  the 
present  work,  the  following  sinusoidal  functions  originally 
proposed  by  Ito  and  Chou  [9]  are  used: 


h  f  .  2n  .  2tt 
Ztop(x,  y)  =  -  I  sm— x  +  sin— y 


(1) 


h 


2tt 


2tt 


Zbot(x,  y)  =  --  sin— -x  +  sin—y 


(2) 


As  pointed  out  earlier,  fiber  tows  have  typically  a 
near-elliptical  cross-section  and  hence  Eqs.  (1)  and  (2)  only 
approximate  the  actual  tow  cross-section  shape.  Neverthe¬ 
less,  they  are  used  in  the  present  work  since  they  greatly 
simplify  permeability  calculations  in  the  distorted  fabric 
and  are  generally  considered  as  a  good  approximation  for 
the  actual  tow  cross-section  shape. 

A  simple  examination  of  Fig.  2  shows  that  within  a 
single-layer  orthogonal  plain-weave  fabric  unit  cell,  one 
can  identify  three  distinct  domains: 


2.7.  Fabric  architecture 

In  this  work,  only  (un- sheared  and  sheared)  balanced 
orthogonal  plain-weave  fabric  is  considered.  Due  to  the 
in-plane  periodicity,  the  fabric  architecture  can  be  repre¬ 
sented  using  a  unit  cell.  The  entire  orthogonal  plain- weave 
fabric  can  then  be  obtained  by  repeating  the  unit  cell  in  the 
in-plane  (x-  and  y-directions).  A  schematic  of  one  quarter 
of  a  plain-weave  unit  cell  with  the  appropriate  denotation 
for  the  system  dimensions  are  shown  in  Fig.  2.  In  a  typical 
plain- weave  fabric,  the  fabric  thickness  ( h )  to  the  quarter  cell 
in-plane  dimension  (L)  ratio,  h/L,  is  small  (0.01-0.1),  while 
the  tow  cross-section  is  nearly  elliptical  in  shape  with  a  large 


•  the  top  channel,  region  T:  ztop  <  z  <  h/ 2; 

•  the  fabric,  region  F:  zbot  <  z  <  ZtoP; 

•  the  bottom  channel,  region  B:  —  h/2  <  z  <  zbot- 

Regions  T  and  B  contain  only  the  resin,  while  region  F 
contains  both  the  fiber  tows  and  the  resin.  The  resin  flow 
through  a  unit  cell  is  analyzed  in  the  present  paper  by  first 
considering  the  flow  within  the  three  regions  separately  and 
then  utilizing  the  matching  boundary  conditions  which  en¬ 
sure  continuity  in  the  pressure  and  the  velocity  components 
across  the  contact  surfaces  of  the  adjacent  regions.  The  resin 
is  considered  as  a  Newtonian  (constant  density)  fluid.  The 
flows  within  the  top  and  the  bottom  channels  are  assumed  to 
be  of  a  creeping  nature  (i.e.  the  inertial  effects  are  neglected) 


Fig.  2.  Schematic  of  one  quarter  of  the  unit  cell  for  a  one-layer  balanced  orthogonal  plain-weave  fabric. 
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while  the  flow  within  the  fabric  is  assumed  to  be  governed 
by  the  Darcy’s  law  (a  velocity  versus  pressure  gradient  re¬ 
lation  which  eliminates  the  need  for  use  of  the  momentum 
conservation  equations). 


where  zUpp(u  y)  andziowCu  y)  are  mathematical  expressions 
for  the  upper  and  the  lower  surfaces  of  the  channels  and  the 
associated  transverse  velocities,  w\Zupp  and  w\zlow,  are  given 
by  the  appropriate  boundary  conditions  discussed  later. 


2.2.  Governing  equations 


2.2.1.  Flow  within  the  top  and  the  bottom  channels 

Under  typical  fabric  infiltration  conditions,  the  resin  flow 
within  the  regions  T  and  B  can  be  considered  as  a  creeping 
flow  in  which  inertial  effects  are  negligibly  small  in  compar¬ 
ison  to  the  viscous  effects.  Under  such  conditions,  at  con¬ 
stant  temperature,  the  resin  flow  can  be  described  by  the 
Stokes  equations  as: 


dp  /  dzu  dzu  d zu 

^ ^  dx2  3y2  dz2  dxdy 


dx 


d2V  d2W 

+ 


dxd  Z 


=  0  (3) 


2.2.2.  Flow  within  the  fiber  tows 

The  resin  flow  through  the  fabric  is  described  in  the 
present  work  using  the  Darcy’s  law  for  an  anisotropic  porous 
medium  as: 


v  = 


1  /  dp  dp  dp 

-lKxxf  +  Kxyf  +  Kxzf 
rj  \  dx  dy  dz 

l  f  dp  dp  dp 

-AKyxyX  +  KyyTy  +  KyzK 


1  /  dp  dp  dp 

-[Kvc^  +  Kzyi^  +  Kzzi- 
rj  \  dx  dy  dz 


(8) 

(9) 

(10) 


dp 

-r-  +  *? 
dy 

dp 

a  +  71 

dz 


d2v 


d2v 


d2u 


2~  ~  -  -  32f 

3y2  dx2  dz2  +  dxdy 


d2w  \ 

H - =  0 

dydzj 


(4) 


(  d2W 
dz2 


d2W 

+  5V  + 


d2W 


d2U 


dy2  dzdx 


+ 


d2v 

dydz 


du  dv  dw 

- 1 - I - =0 

dx  dy  dz 


0  (5) 


(6) 


where  p  is  the  pressure,  u ,  v  and  w  are  respectively  the  x-, 
y-  and  z-components  of  the  resin  velocity  and  q  the  resin 
viscosity. 

Following  the  procedure  of  Simacek  and  Advani  [7], 
which  involves  non-dimensionalization  of  the  governing 
equations,  and  the  use  of  the  conditions:  h/L  1  and 
WL/Uh  ~  1  (U  and  W  are  (mean)  in-plane  and  transverse 
resin  velocity  magnitudes,  respectively),  Eqs.  (3)-(6)  can 


be  simplified  to  yield: 

dp  d2u 

~Vx+V^=Q 

O') 

dp  d2v 

~^  +  n&=0 

(40 

CL5  I 

II 

o 

(50 

du  dv  dw 
—  T  —  T  —  =  0 
dx  dy  dz 

(60 

Eqs.  (3/)-(6/)  are  generally  referred  to  as  “two-dimensional 
lubrication-flow  equations”  in  which  the  pressure  variation 
in  the  z-direction  is  negligibly  small.  However,  in  contrast 
to  the  traditional  lubrication  models,  the  transverse  velocity 
w  (the  velocity  in  the  z-direction)  is  generally  not  zero  (or 
constant)  in  the  present  case  and,  consequently,  the  last  term 
on  the  left  hand  side  of  the  continuity  equation,  Eq.  (60, 
does  not  vanish.  Nevertheless,  this  term  can  be  eliminated 
by  integrating  Eq.  ( 6' )  in  the  z-direction  to  yield: 
fz upp  ( du  dv\ 

L  U+wdz+“,|'~-“’i»-=o  m 


du  dv  dw 
—  T  —  T  —  =  0 
dx  dy  dz 


(11) 


where  Kxx,  Kyy,  Kzz,  Kxy=Kyx,  KXZ=KZX,  and  Kyz=Kzy  are 
the  components  of  the  symmetric  tow  permeability  tensor. 

Eqs.  (8)— (1 1)  can  be  simplified  under  the  following  as¬ 
sumptions:  (a)  Kxx  =  Kyy  =  and  (b)  Kyz  =  Kzx  = 
0.  The  first  assumption  is  not  typically  fully  justified  since 
the  longitudinal  components  of  the  permeability  ( Kxx  and 
Kyy)  are  generally  larger  (up  to  an  order  of  magnitude)  than 
the  transverse  component  ( Kzz )  of  the  permeability.  How¬ 
ever,  this  assumption  greatly  simplifies  the  computational 
procedure  and,  for  simple  fabric  geometries,  it  is  found, 
in  the  present  work,  that  the  results  are  different  by  only 
1-2%  relative  to  their  more  accurately  determined  counter¬ 
parts  corresponding  to  Kxx/Kzz  =  Kyy/Kzz  =  10.  The  sec¬ 
ond  assumption,  on  the  other  hand,  is  generally  expected 
to  be  valid  for  at  least  two  reasons:  (a)  for  the  orthogonal 
plain- weave  architecture  of  the  fabric,  the  material  trans¬ 
verse  principal  direction  is  expected  to  be  essentially  coin¬ 
cident  with  the  global  z-axis;  and  (b)  the  second  assumption 
is  valid  whenever  the  first  assumption  is  valid.  Again  fol¬ 
lowing  the  procedure  of  Simacek  and  Advani  [7],  which  in¬ 
volves  non-dimensionalization  of  the  governing  equations, 
and  the  use  of  the  conditions:  h/L  ^  1  and  WL/Uh  ~  1, 
Eqs.  (8)— (1 1)  become: 


u  =  0  (80 

v  =  0  (90 


Kzz/ )p 

rj  3z 


(100 

(HO 


Eqs.  (80-(ll0  indicate  that  the  only  non-zero  component 
of  the  resin  velocity  within  the  fabric  is  the  one  in  the 
z-direction  and  that,  at  given  values  of  the  in-plane  x-  and 
y-coordinates,  this  component  of  the  velocity  does  not  vary 
in  the  z-direction. 
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2.3.  Boundary  conditions 

The  following  boundary  conditions  are  used  for  the  resin 
flow  problem  in  the  three  regions: 

•  no  slip  (u  =  v  =  w  =  0)  at  the  mold  walls,  z  =  ±(h/2); 

•  at  the  fabric/channels  contact  surfaces,  ztop  and  Zbot>  the 
velocities  and  the  pressure  continuity  are  assumed,  i.e.: 

0(inT)|Ztop=0(inF)|Ztop  (12) 

0(inB)|Zbot  =  0(inF)| 

£bot  (13) 

where  0  =  p,  u,  v,  or  w. 

It  should  also  be  noted  that,  as  established  in  the  previous 
section,  u( inF)  =  i>(inF)  =  0.  In  addition,  the  definition  of 
the  (x  and  y)  in-plane  boundary  conditions  is  deferred  until 
the  final  system  of  equations  is  derived  (the  next  section). 


2.4.  The  final  system  of  equations 


The  resin  velocities  in  the  two  channels  can  be  obtained  by 
integrating  twice  Eqs.  (12)  and  (13),  and  using  the  boundary 
conditions  given  by  Eqs.  (12)  and  (13)  to  determine  the 
integration  constants.  This  procedure  yields: 


dp/dx  (  h 

u  =  -  I  Z  -  -  )  (z  -  "top) 


V 

dp/dy 

t] 

dp/dx 

T] 

dp/dy 
1 1 


(z- 


h  . 

2  I  (z  -  ztop) 


'z+\)(z 


'  "bot) 


/  h 

(Z+2 


(z  £bot) 


in  regional  T 


in  region  B 


(14) 


(15) 


The  subsequent  equations  can  be  simplified  by  introducing 
the  following  expressions:  hT(x,  y)  =  h/2— ztop,  hB(x,  y)  = 
Zbot  -  h/2,  hF(x ,  y)  =  Ztop  -  Zbot?  which  denote  the  height 
fields  of  the  top  channel  and  the  bottom  channels  and  the 
thickness  field  of  the  fabric,  respectively. 

Substitution  of  Eqs.  (14)  and  (15)  into  the  integrated  form 
of  the  continuity  equation,  Eq.  (7),  for  the  two  channels 
yields: 


w  \z=h/2  w  lz=ztop 

_  J_  /8((hT)3(dpJ/dx))  9((/iT)3QpT/9y)) 

6q  \  dx  dy 

w  l^=^bot  —  w  \z=—h/2 

_  J_  /9((/tB)3(9pB/9x))  +  9((/;B)3(9pB/9 y)) 
6rj  \  dx  dy 


=  0 
(16) 

=  0 
(17) 


where  superscripts  T  and  B  are  used  to  denote  the  quantities 
pertaining  to  the  top  and  the  bottom  channels. 


The  first  two  terms  on  the  left  hand  side  of  Eqs.  (16) 
and  (17)  are  defined  by  the  boundary  conditions  discussed 
earlier  as: 


T 

W 
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II 
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W 

1 Z — Zbot 

qhF 

Consequently  Eqs.  (16)  and  (17)  can  be  rewritten  as 

-  P) 

r)hfi 

_  J_  /  9((/iT)3(9pT/9x))  +  9((/iT)3(9j?T/9y)) 

6  x]  \  dx  dy 

(22) 

^(FB  ~  P) 

_  J_  / 9((/iB)3(9pB/9x))  9((/iB)3(9pB/9y)) 

6r]  \  dx  dy 

(23) 

Eqs.  (22)  and  (23)  represent  the  final  system  of  equations 
consisting  of  two  coupled  linear  elliptic  partial  differential 
equations  with  the  pressures  pT  and  pB  as  dependent  vari¬ 
ables.  To  solve  these  equations,  boundary  conditions  along 
the  (x-y)  in-plane  boundaries  of  the  unit  cell  must  be  pre¬ 
scribed.  For  the  un- sheared  balanced  plain- weave  fabric  ar¬ 
chitecture  in  which  the  unit  cell  boundaries  are  the  lines  of 
geometrical  symmetry,  a  fixed  pressure  gradient  can  be  en¬ 
forced  in  one  principal  direction  while  requiring  periodicity 
in  the  pressure  distribution  in  the  direction  normal  to  the 
direction  in  which  the  pressure  gradient  is  prescribed.  This 
type  of  boundary  conditions  is  generally  used  since  it  en¬ 
ables  determination  of  the  off-diagonal  (/C^,  Kyz  and  Kxz) 
components  of  the  effective  preform  permeability  tensor.  A 
more  detailed  discussion  of  the  in-plane  boundary  conditions 
is  given  later  in  the  context  of  the  effect  of  fabric  shearing 
on  the  choice  of  in-plane  boundary  conditions. 

The  final  system  of  partial  differential  equations,  Eqs.  (22) 
and  (23),  contains  the  thickness  fields:  hT(x,  y),  hB(x,  y) 
and  hF(x,  y).  These  fields  are  defined  in  the  present  work 
using  the  analytical  expressions  for  the  top  and  the  bottom 
surfaces  of  the  fabric  preform  (Eqs.  (1)  and  (2)).  These  ex¬ 
pressions  are  generally  considered  as  reasonably  good  ap¬ 
proximations  of  the  actual  orthogonal  plain- weave  fabric  ar¬ 
chitecture  with  a  near  elliptical  cross-section  area.  It  should 
be  noted,  however,  that  over- simplification  of  the  fabric  ar¬ 
chitecture  (e.g.  using  square  or  circularly  shaped  tows)  may 
lead  to  erroneous  results  and  must  be  avoided.  In  general, 
the  thickness  fields  can  be  constructed  using  direct  exper¬ 
imental  measurements  such  as  quantitative  metallographic 
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analysis  of  consolidated  and  sectioned  parts,  e.g.  [11],  and 
through  the  use  of  computerized  image  analysis  of  the  fab¬ 
ric  surface,  e.g.  [12].  The  second  of  these  two  methods  is 
quite  appealing  since  the  image  conversion  procedure  can 
be  directly  coupled  with  the  solution  scheme  for  Eqs.  (22) 
and  (23). 

Due  to  complexity  in  the  hT(x,  y),  hB(x,  y )  and  hF (x,  y ) 
functions,  Eqs.  (22)  and  (23),  cannot  be  solved  analytically. 
However,  finding  the  numerical  solution  to  Eqs.  (22)  and 
(23)  is  relatively  straightforward.  In  the  present  work,  MAT- 
LAB  general-purpose  mathematical  package  [10]  and  a  fi¬ 
nite  difference  method  are  used  to  solve  Eqs.  (22)  and  (23). 

Once  Eqs.  (22)  and  (23)  are  solved,  the  resulting  pressure 
fields  can  be  used,  in  conjunction  with  Eqs.  (14)  and  (15),  to 
compute  the  corresponding  in-plane  velocity  fields  in  the  two 
channels.  Integration  of  these  velocity  fields  over  the  side 
boundaries  of  the  quarter  unit  cell  then  enables  determine  of 
the  total  resin  flow  rate,  Q  =  [Qx,  Qy],  through  the  quarter 
unit  cell  in  the  two  principal  direction.  The  components  of 
the  effective  in-plane  preform  permeability,  and 

Kfy  are  then  computed  using  the  two-dimensional  Darcy’s 
law  and  the  known  imposed  values  of  the  pressure  gradient. 


2.5.  Application  of  the  model  to  the  multi-layer  fabric 


The  model  developed  thus  far  pertains  to  a  single-layer 
fabric  preform.  In  typical  RTM  and  VARTM  processes, 
the  preforms  may  contain  several  fabric  layers.  In  such 
multi-layer  preforms,  nesting  and  compaction  generally 
have  a  significant  effect  and  must  be  included  when  pre¬ 
dicting  preform  permeability.  Numerous  experiments,  e.g. 
[13,14],  confirmed  that  permeability  varies  with  a  number 
of  layers. 

The  single-layer  model  developed  in  the  previous  sec¬ 
tion  can  be  readily  extended  to  a  multi-layer  preform.  A 
schematic  of  two  types  of  two-layer  plain- weave  fabric  pre¬ 
forms  is  given  in  Fig.  3(a)  and  (b).  The  two  types  are  gen¬ 
erally  referred  to  as  “in-phase”  and  “out-of-phase”  fabric 
architectures  or  laminates.  In  the  case  of  an  ft-layer  fabric 
preform,  if  the  channels  are  labeled  using  consecutive  in¬ 
tegers  (with  the  bottom  channel  being  denoted  as  channel 
“1”),  the  analytical  procedure  for  a  single-layer  fabric  pre¬ 
form  used  in  the  previous  section  yields  (n  +  1)  coupled 
elliptical  partial  differential  equations  with  n  +  1  unknown 
pressures  p^\  p^2\  . . .  ,  p(nFl>}  as: 

-%r(P(2)  ~  PW)  ~  fv((hm)3Vpm)  =  0  (24) 

77/r 1  or] 


5^{p(‘)  _  /+D)  _  (pO'+D  -  p(0) 


r]hFi- 


6r] 


r\hFi 

V((ft(0)3Vp(0)  =  0,  i  =  2,3,...  ,n 


(25) 


K 


4r(P(n+l)  ~  P(n))  ~  —V((h{n+l))Vpin+l))  =  0  (26) 

r]htn  677 


► 


X 


Fig.  3.  x-z  section  of  a  quarter  of  the  unit  cell  for  (a)  an  in-phase  and 
(b)  an  out-of-phase  two-layer  orthogonal  plain-weave  fabric. 


where  hFi  (i  =  1,  . . .  ,n )  denotes  the  thicknesses  of  the  ith 
fabric  layer  (numbered  starting  from  the  bottom  of  the  mold) 
and  h®  (i  =  1, . . .  ,  n  +  1)  are  the  heights  of  the  inter-fabric 
or  tool/fabric  resin  channels  (also  numbered  starting  from 
the  bottom  of  the  mold).  The  system  of  equations  defined 
by  Eqs.  (24)-(26)  is  solved  using  the  same  computational 
procedure  used  for  the  one-layer  fabric  preform. 

2.6.  Shear-induced  fiber  volume  fraction  correction  for 
permeability 

When  the  fabric  is  sheared  in  the  v-direction,  as  shown 
in  Fig.  7(b),  weft  tows  are  rotated  but  remain  stress  free. 
Consequently,  the  dimension  of  the  fabric-preform  unit  cell 
in  the  y-direction  is  altered  causing  a  change  in  the  effec¬ 
tive  fiber  volume  fraction  in  the  unit  cell.  This  change  in  the 
fiber  volume  fraction  can  have  a  significant  effect  on  pre¬ 
form  permeability  at  large  shear  angles  and,  hence,  must  be 
taken  into  account.  The  procedure  described  below  is  used  to 
correct  permeabilities  in  the  sheared  fabrics  obtained  using 
the  original  model  of  Simacek  and  Advani  [7]  as  reviewed 
in  Sections  2.4  and  2.5. 
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Fig.  4.  (a)  Resin  channels  height  and  (b)  fabric  thickness  fields  in  an  un-sheared  one-layer  orthogonal  plain-weave  fabric  preform. 


To  quantify  the  permeability  correction  described  above, 
the  Kozeny-Carman  relation,  e.g.  [20],  for  permeability  of 
the  porous  media  with  a  fibrous  architecture  is  used.  Ac¬ 
cording  to  this  relation,  permeability  of  such  media  is  given 
by: 


K  = 


F(i-/)3 

cfi 


(27) 


where  r  and  /  are  the  fiber  radius  and  the  fiber  vol¬ 
ume  fraction,  respectively,  while  c  is  a  fibrous -medium 
architecture-dependent  constant. 

When  the  fabric  preform  is  sheared  in  the  v-direction  by 
an  angle  0 ,  the  fiber  volume  fraction  in  fabric  tows  changes 
as: 


fe 


fo 

sin  (90  —  0) 


(28) 


where  the  angle  0  is  given  in  degrees  and  the  subscripts  0 
and  6  are  used  to  denote  the  value  of  a  respective  quantity  in 
the  un-sheared  fabric  and  in  the  fabric  sheared  by  an  angle 
6,  respectively. 

To  account  for  a  shear-induced  change  in  the  fiber  volume 
fraction,  the  permeability  values  for  sheared  fabric  preforms 
obtained  using  the  models  described  in  Sections  2.4  and  2.5, 
should  be  multiplied  by  the  following  correction  factor: 


/o2(l  -  fo)3 


(1  -  /o)3/*2 


(29) 


3.  Results  and  discussion 

3.1.  Un-sheared  single-layer  plain-weave  fabric  preforms 

The  model  developed  in  Section  2.4  is  used  in  this  section 
to  analyze  the  pressure  distribution  within  the  un-sheared 
single-layer  balanced  orthogonal  plain-weave  quarter  unit 


cell.  Due  to  the  symmetry  of  the  unit  cell  with  respect  to  the 
z  =  0  plane,  the  pressure  distributions  within  the  top  and 
the  bottom  resin  channels  are  identical  and,  hence,  there  is 
no  transverse  flow  of  the  resin  through  the  fabric  preform. 
Also,  as  established  in  Section  2.4,  there  is  no  variation  of 
the  pressure  in  the  z-direction  within  the  channels.  The  vari¬ 
ation  of  the  top-  and  bottom-channel  heights  and  of  the  fab¬ 
ric  thickness  in  the  x-y  plane  within  a  quarter  unit  cell,  used 
as  input  in  the  present  analysis,  are  shown  in  Fig.  4(a)  and 
(b),  respectively.  The  variation  of  the  pressure  in  the  x-y 
plane  within  the  resin  channels  of  a  quarter  unit  cell  for  the 
fixed  pressure  drop  of  1.0  x  105  in  the  v-direction  is  shown 
in  Fig.  5.  The  pressure  distribution  (or  more  precisely  its 


Distance  in  x-direction,  m 


Fig.  5.  Pressure  distribution  in  the  x-y  plane  within  a  resin  channel  in 
the  case  of  an  un-sheared  single-layer  balanced  orthogonal  plain-weave 
fabric  preform. 
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gradient)  at  a  given  x-y  location  correlates  inversely  with  the 
local  height  of  the  resin  channel  in  order  to  satisfy  the  con¬ 
tinuity  equation.  It  should  be  also  noted  that  due  to  the  sym¬ 
metry  of  the  fabric  geometry  with  respect  to  the  quarter  unit 
cell  boundaries  normal  to  the  y-direction,  zero-flux  (i.e.  zero 
pressure  gradient)  conditions  are  found  in  the  y-direction  at 
these  boundaries. 


3.2.  Effect  of  the  number  of  layers  in  un- sheared 
plain-weave  fabric  preforms 

The  model  developed  in  Section  2.5  is  used  in  the  present 
section  to  predict  permeability  of  the  balanced  un-sheared 
single-  and  multi-layer  orthogonal  plain-weave  fabric  archi¬ 
tectures.  In  all  the  calculations  carried  out  in  this  section,  as 
well  as  in  the  calculations  carried  out  in  the  previous  sec¬ 
tion,  the  following  unit  cell  parameter  and  one-layer  fabric 
thickness  values  are  used:  L\  =  L2  =  L  =  0.01m  and 
h  =  0.001  m.  Also  the  transverse  permeability  of  the  fiber 
tows  is  set  to  a  typical  (fixed)  value,  Kzz  =  1  x  10-10  m2. 

To  determine  the  effect  of  the  number  of  fabric  layers 
on  the  effective  permeability,  the  model  developed  in  the 
previous  section  is  used  for  the  cases  of  1-,  2-,  3-,  5-,  10- 
and  20-layer  in-phase  orthogonal  balanced  plain- weave  fab¬ 
ric  preforms  in  the  absence  of  layer  nesting.  The  results 
of  this  calculation  are  presented  in  Fig.  6.  These  results 
show  that  as  the  number  of  layers  increases,  the  permeabil¬ 
ity  rises  but  at  an  ever  decreasing  rate  so  that  in  fabric  pre¬ 
forms  with  10  or  more  layers,  the  effect  of  the  number  of 
layers  on  permeability  becomes  insignificant.  This  finding 
can  be  easily  rationalized  by  recognizing  that  as  the  num¬ 
ber  of  layers  in  the  fabric  increases,  the  effect  of  the  bottom 
and  the  top  resin  channels  which  are  more  restrictive  to  the 
fluid  flow  (and  thus  reduce  effective  preform  permeability) 
decreases. 


Fig.  6.  The  effect  of  the  number  of  fabric  layers  on  the  effective  perme¬ 
ability  of  an  un-sheared  balanced  plain-weave  fabric  preform. 


3.3.  Effect  of  fabric  shear  on  permeability 

As  pointed  out  earlier,  when  the  fabric  preform  is  forced 
to  conform  to  the  ridges  and  recesses  of  a  mold,  it  may 
locally  undergo  shear  deformation.  Such  deformation  can 
significantly  affect  local  permeability  of  the  preform.  As 
shown  in  Fig.  7,  when  a  balanced  square-cell  plain- weave 
fabric  is  sheared,  two  important  factors  must  be  consid¬ 
ered:  (a)  the  unit  cell  size  increases  and  to  make  the  calcu¬ 
lations  of  preform  permeability  manageable,  the  shear  an¬ 
gle  a  =  tan ~l(m/n)  is  generally  allowed  to  take  only  the 
values  corresponding  to  relatively  small  integers  m  and  n\ 
and  (b)  the  boundaries  of  the  unit  cell,  unlike  the  case  of 
the  initial  square-shape  unit  cell,  are  no  longer  the  lines  of 
symmetry  of  the  fabric  structure.  Consequently,  the  bound¬ 
ary  conditions  imposed  along  the  boundaries  of  the  unit  cell 
have  to  be  modified  relative  to  those  used  in  the  case  of  the 
un-sheared  unit  cell.  For  instance,  if  a  fixed  pressure  drop 
is  applied  in  the  v-direction,  the  symmetry  conditions  along 
the  unit  cell  boundaries  normal  to  the  y-coordinate  require 
that  zero  pressure-gradient  boundary  conditions  be  applied 
in  the  y-direction.  In  the  case  of  a  sheared  fabric  preform, 
on  the  other  hand,  the  unit  cell  boundaries  normal  to  the 
y-direction  are  not  any  longer  the  lines  of  symmetry  of  the 
fabric  architecture  and,  hence,  only  the  periodic  boundary 
condition  (the  corresponding  pressure  values  along  the  two 


Fig.  7.  Effect  of  fabric  shearing  on  the  size  of  the  quarter  unit  cell  (denoted 
using  heavy  dashed  lines)  in  balanced  plain- weave  fabric  architectures: 
(a)  un-sheared  fabric;  (b)  fabric  sheared  by  an  angle  a  =  tan-1  (1/3)  in 
the  x-direction. 
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Fig.  8.  The  effect  of  shear  on  permeability  of  a  single-layer  balanced 
orthogonal  plain-weave  fabric  preform. 

unit-cell  boundaries  normal  to  the  y-direction  are  identical) 
can  be  applied. 

The  effect  of  shear  deformation  (measured  by  the  mag¬ 
nitude  of  the  shear  angle  a)  on  the  effective  permeabil¬ 
ity  of  single-layer  plain-weave  fabric  preforms  is  displayed 
in  Fig.  8.  An  example  of  the  variation  of  the  top-  and 
bottom-channel  heights  and  of  the  fabric  thickness  in  the 
x-y  plane  within  a  quarter  unit  cell,  used  as  input  in  the 
present  analysis,  are  shown  in  Figs.  9(a)  and  (b),  respec¬ 
tively.  For  comparison,  the  experimental  values  of  preform 
permeabilities  obtained  in  Ref.  [8]  are  also  shown  in  Fig.  8. 
While  the  agreement  between  the  corresponding  computa¬ 
tional  and  the  experimental  values  is  only  fair,  the  effect  of 
shear  deformation  on  preform  permeability  appears  to  be 
quite  well  predicted  by  the  model.  In  addition,  the  corre¬ 
sponding  computed  values  of  the  in-plane  off-diagonal  (Kxy 
and  Kyx)  elements  of  the  effective  permeability  are  very 
close  as  required  by  symmetry  of  the  orthogonal  plain-weave 
fabric  architecture. 


Fig.  10.  A  typical  compaction-pressure  vs.  preform-thickness  curve  for  a 
plain-weave  fabric  architecture. 

3.4.  Effect  of  preform  compaction  on  permeability 

When  the  fabric  is  subjected  to  compression  during  mold 
closing  in  the  RTM  process  or  during  evacuation  of  the  vac¬ 
uum  bag  in  the  VARTM  process,  it  undergoes  a  number  of 
changes  such  as:  the  cross-section  of  the  fiber  tow  flattens, 
the  pores  and  gaps  between  the  fibers  inside  tows  as  well  as 
between  individual  tows  are  reduced,  the  tows  undergo  elas¬ 
tic  deformation,  inter-layers  shifting  (nesting),  etc.  A  typical 
compression-pressure  versus  preform  thickness  curve  for  a 
woven  fabric  is  depicted  in  Fig.  10  [15].  The  curve  shown  in 
Fig.  10  has  three  distinct  parts:  two  linear  and  one  nonlinear. 
In  the  low-pressure  linear  and  the  nonlinear  portions  of  the 
pressure  versus  thickness  curve,  preform  compaction  is  dom¬ 
inated  by  a  reduction  of  the  pore  and  the  gap  sizes  between 
the  fibers  in  tows.  In  the  high-pressure  linear  region  of  the 
pressure  versus  thickness  curve,  on  the  other  hand,  preform 


Fig.  9.  (a)  Resin  channels  height  and  (b)  fabric  thickness  fields  in  a  one-layer  orthogonal  plain- weave  fabric  preform  subjected  to  shear  in  the  x-direction 
by  an  angle  of  a  =  tan-1  (1/3). 
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compaction  involves  mainly  tow  bending  and  nesting.  Typ¬ 
ical  liquid  molding  processes  such  as  RTM  or  VARTM  in¬ 
volve  pressures  which  correspond  to  the  high-pressure  lin¬ 
ear  pressures  versus  thickness  region.  Hence,  the  effect  of 
fabric  compaction  on  permeability  of  the  fabric  preform  as¬ 
sociated  with  the  high-pressure  linear  compaction  regime  is 
investigated  in  this  section. 

To  quantify  the  effect  of  preform  compaction  (in  the 
high-pressure  linear  region)  on  permeability  of  the  balanced 
orthogonal  plain- weave  fabric,  the  beam-bending-based 
micro-mechanical  model  developed  in  a  series  of  papers  by 
Chen  and  Chou  [16-18]  is  utilized  in  the  present  work.  The 
model  of  Chen  and  Chou  [16-18]  is  based  on  a  number  of 
well-justified  assumptions  such  as:  (a)  the  fabric  is  consid¬ 
ered  to  extend  indefinitely  in  the  x-y  plane  and,  hence,  can 
be  represented  using  the  unit  cells  such  as  the  one  shown 
in  Fig.  2;  (b)  tows  in  the  fabric  are  treated  as  a  transversely 
isotropic  solid  material;  (c)  the  fabric  is  subjected  to  the 
compaction  pressure  only  in  the  through-the-thickness  di¬ 
rection,  and  can  freely  adjust  its  shape  in  the  x-y  plane; 
(d)  since  the  compaction  analyzed  corresponds  to  the  high- 
pressure  linear  region,  no  voids  or  gaps  are  assumed  to 
exist  between  the  fibers  in  tows  or  between  the  tows;  (e) 
during  fabric  compaction,  the  cross-section  area  of  the 
tows  is  assumed  to  remain  unchanged  but  the  shape  of  the 
cross-section  undergoes  a  change;  and  (f)  as  compaction 
proceeds,  the  deformation  of  the  tows  leads  to  an  increase 
in  the  effective  volume  fraction  of  the  fibers  in  the  fabric 
and,  in  the  limit  of  complete  compaction  of  the  tows,  the 
volume  fraction  of  the  fibers  in  the  fabric  becomes  equal  to 
that  in  the  individual  tows. 

In  order  to  derive  a  relationship  between  the  reduction 
in  the  fabric  thickness,  the  effective  volume  fraction  of  the 
fibers  and  various  distributions  and  magnitudes  of  the  ap¬ 
plied  compaction  pressure,  Chen  and  Chou  [16-18]  applied 
a  simple  procedure  from  the  solid  mechanics  beam  theory. 
Toward  that  end,  the  one-quarter  unit  cell  shown  in  Fig.  2  is 
first  simplified  by  replacing  the  two  warp  and  the  two  weft 
tows  with  four  beams.  Next  based  on  the  symmetry  of  the 
simplified  model,  it  is  shown  that  the  problem  can  be  further 
simplified  using  a  single  beam  and  the  appropriate  distribu¬ 
tion  of  the  applied  and  contacting  pressures  (Fig.  11).  The 
model  of  Chen  and  Chou  [16-18]  is  utilized  in  the  present 
work  to  compute  the  effect  of  the  compaction  pressure  on 
the  channel  heights  ( hT(x ,  y)  and  hB(x,  y))  and  on  the  fabric 
thickness,  hF (x,  y )  fields.  These  fields  are,  in  turn,  used  in 
the  lubrication  model  presented  in  the  Section  2.4  to  quan¬ 
tify  the  effect  of  fabric  compaction  on  the  effective  perme¬ 
ability  of  a  one-layer  orthogonal  plain- weave  fabric. 

The  effect  of  the  compaction  force  applied  to  the  upper 
and  the  lower  (rigid  and  flat)  molds  on  the  effective  perme¬ 
ability  of  a  one-layer  orthogonal  plain-weave  fabric  is  dis¬ 
played  in  Fig.  12.  In  these  calculations,  the  Young’s  modulus 
is  assigned  a  value  of  22  GPa  and  a  sinusoidal  distribution  of 
the  applied  and  the  contacting  pressures  is  assumed  [18].  An 
example  of  the  variation  of  the  top-  and  the  bottom-channel 


Applied  Pressure 


Distribution 


Fig.  11.  A  schematic  of  the  pressure  distribution  on  a  curved  beam  used 
in  the  calculation  of  permeability  of  un-sheared  one-layer  orthogonal 
plain- weave  fabrics. 

heights  and  of  the  fabric  thickness  in  the  x-y  plane  within 
a  quarter  unit  cell,  used  as  input  in  the  present  analysis,  are 
shown  in  Fig.  13(a)  and  (b),  respectively.  For  comparison, 
the  experimental  results  reported  by  Sozer  et  al.  [19]  are  also 
shown  in  Fig.  12.  It  is  seen  that  a  reasonably  good  agree¬ 
ment  exists  between  the  both  the  magnitude  of  the  predicted 
preform  permeability  and  its  change  with  the  applied  com¬ 
paction  force. 

3.5.  Effect  of  layer  nesting 

As  mentioned  earlier  shifting  of  fabric  layers  followed 
by  their  more  compact  packing  (the  phenomenon  generally 
referred  to  as  layers  “nesting”)  can  have  a  major  effect  on 
the  effective  fiber  density  in  the  preform  and,  hence,  on 


Fig.  12.  Effect  of  compaction  (represented  by  the  magnitude  of  the 
compaction  force)  on  permeability  of  a  one-layer  un-sheared  orthogonal 
plain-weave  fabric  preform. 
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Fig.  13.  (a)  Resin  channels  height  and  (b)  fabric  thickness  fields  in  an  un-sheared  one-layer  orthogonal  plain-weave  fabric  preform  subject  to  a  total 
compressive  force  of  10.3  N  via  rigid,  flat  upper  and  lower  tool  surfaces. 


permeability  of  the  preform.  Nesting  of  the  fabric  layers  can 
particularly  take  place  under  high  applied  pressures  which 
are  sufficient  to  overcome  inter-tow  friction.  The  thickness 
reduction  in  balanced  orthogonal  plain-weave  fabrics  whose 
geometry  is  represented  by  Eqs.  (1)  and  (2),  due  to  layers 
nesting,  has  been  analyzed  by  Ito  and  Chou  [9]  who  de¬ 
rived  the  following  relation  for  the  fabric  thickness  reduction 
caused  by  nesting: 


and  fabric  permeability  at  4>x  =  0-y  =  0.  As  expected,  fabric 
nesting  gives  rise  to  the  reduction  in  fabric  permeability. 
Furthermore,  for  the  case  of  a  out-of-phase  laminate  fabric 
(4>x  =  (py  =  ±(7t/2)),  fabric  permeability  is  only  about  30% 
of  its  value  in  the  in-phase  laminate  fabric.  This  finding  is 
in  excellent  agreement  with  its  experimentally  counterpart 
reported  by  Sozer  et  al.  [19]. 


A/z  nesting  =  h 


n.  0*  0y 
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IX  7 X 
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10*1  <  |  <  I0y I  <  7X, 

=  <  10*1  <  7t,  =  <  |0y  I  <  7t, 


(30) 


where  (j)x  =  ( 2tx/L)sx  and  (j)y  =  (2jx/L)sy  are  dimensionless 
while  and  sy  are  the  dimensional  relative  shifts  of  the 
adjacent  layers  in  the  x-  and  y-directions,  respectively. 

Two  non-nesting  cases  associated  with  zero  nesting,  re¬ 
duction  in  the  fabric  thickness  can  be  identified:  (a)  0*  = 
(j)y  =  0  which  corresponds  to  the  iso-phase  laminate  case  and 
(b)  cj)x  =  0-y  =  =b(7r/2)  corresponding  to  the  out-of-phase 
laminate  case. 

The  relations  given  in  Eq.  (30)  are  used  in  the  present 
work  to  examine  the  effect  of  layers  nesting  on  fabric  per¬ 
meability.  While,  in  general,  fabric  compaction  during  the 
high-pressure  linear  compaction  stage  can  involve  both  elas¬ 
tic  distortions  (tow  bending)  and  layers  nesting,  the  two 
modes  of  fabric  compaction  are  generally  considered  as  de¬ 
coupled  and  can  be  considered  separately. 

The  effect  of  nesting  (quantified  by  the  magnitudes  of  the 
dimensionless  layer  shifts  in  the  x-  and  the  y-directions,  4>x 
and  0-y,  respectively)  in  a  two-layer  orthogonal  plain- weave 
fabric  is  shown  in  Fig.  14.  The  values  displayed  pertain  to  the 
ratio  of  fabric  permeability  at  the  given  values  of  4>x  and  4>y 


Dimensionless  Shift  in  the  x-Direction 


Fig.  14.  The  effect  of  nesting  on  the  ratio  of  fabric  permeability  at 
the  given  values  of  layer  shifts  in  the  *-  and  y-directions  and  fabric 
permeability  of  an  un-nested  in-phase  laminate  fabric. 
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4.  Conclusions 

Based  on  the  results  obtained  in  the  present  work,  the 
following  main  conclusions  can  be  drawn: 

1.  Effective  permeability  of  the  orthogonal  plain- weave 
fabric  preforms  can  be  determined  computationally 
by  combining  a  lubrication  model  for  the  resin  flow 
through  tool-surface/fabric-tow  and  tow/tow  channels 
with  the  Darcy’s  law  for  the  resin  flow  through  the 
fabric  tows. 

2.  The  computational  approach  presented  in  this  work  en¬ 
ables  assessment  of  the  contribution  that  various  phe¬ 
nomena  such  as  the  mold  walls,  fabric  shearing,  in¬ 
terlayer  shifting  and  restacking  as  well  as  fabric  com¬ 
paction  due  to  the  infiltration  pressure  make  to  orthog¬ 
onal  plain- weave  fabric-preform  permeability. 

3.  While  no  comprehensive  set  of  experimental  data  is 
available  to  fully  test  validity  of  the  present  model, 
the  agreement  of  the  model  predictions  with  selected 
experimental  results  can  be  generally  qualified  as 
reasonable. 
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